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Abstract. Context Population II pulsating variables play a relevant role both as distance indicators and as tracers 
of the properties of old stellar populations. 

Aims In this paper we present an updated and homogeneous pulsational scenario for a wide range of stellar 
parameters typical of BL Her stars i.e., Population II Cepheids with periods shorter than 8 days. Methods To this 
purpose, we adopt a nonlinear convective hydrodinamical code evaluate the stability and full amplitude behaviour 
of an extensive set of BL Her pulsation models. Various assumptions of mass, luminosity and metallicity consistent 
with the most recent evolutionary prescriptions, are adopted. Results We obtain the theoretical instability strip 
for both fundamental and first overtone pulsators and present a detailed atlas of light /radial- velocity curves. 
Some relations for the boundaries of the instability strip and for the dependence of the absolute magnitude on 
period, mass and color, as well as the fundamental period-amplitude relations are derived. Finally, we provide the 
theoretical period-radius relation for BL Her and find that it is in excellent agreement with the empirical relation 
by Burki & Meylan and consistent with the one holding for shorter periods for RR Lyrae stars. 



1. Introduction 

The BL Her variables form a small but interesting group 
of radially pulsating Population II stars with period P in 
the range 0.8-3 days. They are observed in globular clus- 
ters with few RR Lyrae stars and blue horizontal branch 
(HB) morphology and they are brighter than RR Lyrae 
but fainter than Anomalous Cepheids, for a fixed period. 
They are currently interpreted as central He-burning low 
mass (M < 0.6M Q ) stars which, after their main core 
He-burning phase spent in the blue side of the Zero Age 
Horizontal Branch (ZAHB), evolve toward the Asymptotic 
Giant Branch (AGB) crossing the instability strip at lu- 
minosities higher than the RR Lyrae level (see Gingold 
1985; Bono et al. 1997a, Caputo 1998). 
The properties of these objects have recently been re- 
viewed by Wallcrstcin (2002) and Sandage & Tammann 
(2006). From the theoretical point of view, Buchler & 
Moskalik (1992) and Buchler & Buchler (1994) presented 
a pulsational study of models with P < 3 days based on 
a linear and nonlinear radiative analysis. These authors 
found evidence for a resonance between the fundamental 
mode and the second overtone, similar to the one often 
invoked as the explanation of the Hertzsprung progression 
phenomenon in Classical Cepheids (see Bono, Marconi, 
Stcllingwcrf 2000 for a detailed discussion). Buchler & 
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Buchler (1994) also evaluated the pulsation properties of 
first overtone (FO) models and found that the first over- 
tone blue edge (FOBE) was very close (< 100 K) to the 
fundamental (F) one, producing a very narrow regime of 
FO-only pulsation. Moreover, FO models were found to 
show very low pulsation amplitudes but it seemed to be 
hard to discriminate the pulsation mode on the basis of 
their Fourier coefficients. By relying on a nonlinear con- 
vective pulsation analysis, Bono ct al. (1997a) found a 
good agreement between the predicted and the observed 
BL Her instability strip boundaries and suggested that 
these pulsators are pulsating in the F mode with a typical 
mass of 0.52-0.59 M . Moreover these authors predicted 
a Period-Luminosity- Amplitude (PLA) relation in the B 
band. 

However the Bono et al. (1997a) nonlinear convective ap- 
proach was limited to a quite restricted range of stellar 
parameters and adopted an old input physics both in the 
evolutionary and the pulsational computations (see Bono 
& Stellingwerf 1994 and Bono et al. 1997 for details). 
In this paper we present an updated and homogeneous 
pulsational scenario for a wide range of stellar parame- 
ters typical of BL Her stars. In an accompanying paper 
(Di Criscienzo et al. 2006) we will discuss the connection 
between the pulsational scenario and the predictions of 
updated evolutionary models in order to compare theory 
with observations. 
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Table 1. Input parameters of the computed BL Her mod- 
els. An helium abundance Y=0.24 has been adopted. 

Z M/M Q LogL/L FOBE FBE FORE FRE 



0.0001 



0.60 



0.65 



1.95 
2.05 
2.15 
1.91 
2.01 
2.11 



6850 
6750 
6750 
6950 6850 
6750 6850 
6750 



0.001 



0.50 
0.55 

0.65 



2.11 
2.41 
1.81 
1.91 
2.01 
1.81 
1.91 
2.01 



6650 
6350 
6875 6850 
6850 
6750 
7050 6750 
6850 6750 
6650 6850 



0.004 



0.55 



1.81 
1.91 
2.01 



The organization of the present paper is the following: 
in Sect. 2 we present the new computed models. In Sect. 
3 we deal with predicted light/ velocity curves and visual 
amplitudes. Same important relations are given in Section 
4, whereas the conclusions close the paper. 

2. The stellar pulsation models 

To study BL Her stars, we have adopted the same pul- 
sation code and the same physical and numerical as- 
sumptions used for RR Lyrae (Bono & Stellingwcrf 1994, 
Bono, Castellani & Marconi 2000, Marconi et al.2003, Di 
Criscienzo et al. 2004) and Anomalous Cepheid models 
(Marconi et al. 2006). Model sequences have been com- 
puted as one parameter families with constant luminos- 
ity, mass and chemical composition, by varying the effec- 
tive temperature with a step of 100 K. The various as- 
sumptions of mass, luminosity and metallicity are consis- 
tent with the most recent evolutionary prescriptions (see 
Pietrinferni et al. 2004, 2006) and are reported in the first 
three columns of Table 1. 

Linear regression through the model's parameters 
allows us to derive analytical relations connecting the 
period of models to the intrinsic stellar parameters, 
namely mass, luminosity and effective temperature, i.e: 



strip, which depends on the intrinsic parameters of pulsators. 
For each given metallicity, mass and luminosity, we have cal- 
culated the maximum and minimum effective temperature for 
The onset of either fundamental or first-overtone pulsation, as 
5750 reported in the last four columns of Table 1, where FOBE 
5650 (FORE) and FBE (FRE) are the blue (red) limit for first over- 
5550 tone and fundamental pulsator respectively, ft is of importance 
6100 5750 t° note that our calculations confirm the earlier suggestion 
6300 5750 W Bono et al. 1997(a,b,c) that, for a given helium content 
5550 ana - mass, there exists a "intersection" luminosity Lip where 
'_ 5450 FOBE=FBE and that above this luminosity only the funda- 
5150 men t a l mode has a stable nonlinear limit cycle. As shown by 
6400 5650 Criscienzo et al. 2006, the evolutionary models show that 
5550 BL Her variables occur at luminosity higher than this inter- 
5450 section point: therefore, in the following we will consider only 
6650 5750 fundamental models. As for the red boundary of the funda- 
6150 5650 mental mode, it is worth mentioning that it is expected to de- 
6350 5650 P en< i on the efficiency of convection and, in particular, on the 
gggg ~ 5750 vame °f the mixing length parameter I / H p we assume to close 
gggo _ 5650 ^ ne nonlinear system of dynamical and convective equations 
675 q _ 545Q (see Stellingwerf 1982 and Bono & Stellingwerf 1992, 1994 for 
details on the treatment of convection in the adopted hydro- 
dynamical code). To test this occurrence, we have computed 
additional models with l/H p increased from 1.5 to 2.0 and we 
confirm the general trend shown by RR Lyrae (Di Criscienzo et 
al. 2004) and Classical Cepheid models (Fiorentino et al. 2006, 
in preparation) in that the FRE moves toward higher effective 
temperatures as l/H p increases. In addition, we find that this 
effect is less important with increasing the luminosity. In par- 
ticular, for M = O.5OAf and Z=0.001 the FRE blueshift is of 
about 200 K at logLj 'Lq = 2.01 but reduces to less than 100 
K for logL/LQ — 2.41. This occurrence is due to the reduced 
density values, and in turn to the increase of the required su- 
peradiabaticity, in the outer layers of stellar structures with 
low mass and very high luminosity. 



logP = 11.579(±0.015) + 0.893 log ^- - 0.89 log — 

Lq Mq 

- 3.54 log T cff 

for fundamental pulsators and: 

logP = 10.784(±0.003) + 0.806 log ^- - 0.66 log ^~ 

Lq Mq 

- 3.31 log T cff 



(1) 



(2) 



for FO ones. 

It is well known that radial pulsation occurs in a quite 
well-defined region of the HR diagram, the so called instability 



3. Light curves and pulsation amplitudes 

For each investigated model, the non linear analysis provides 
the variation of relevant parameters, namely luminosity, ra- 
dius, radial velocity, effective temperature and gravity, along 
a pulsational cycle. A subsample of bolometric light curves 
(Z=0.001 and M = 0.65M© and the labeled luminosity levels) 
is shown in Figure [T] while Figure [5] exhibits the correspond- 
ing radial velocity curves as a function of pulsation phase. All 
these curves show a large variety of shapes, which is perhaps 
the most striking feature of BL Her models as already sug- 
gested by Moskalik and Buchler (1994). In all the sequences 
(in particular for radial velocity curves), the presence and the 
progression of the main bump is evident in analogy of what 
observed for Classical Cepheids with periods around 10 days 
(the well known Hertzsprung progression). A detailed investi- 
gation of this issue, through the study of Fourier parameters, 
will be addressed in a future paper. 

In order to compare theoretical results with observations, 
the bolometric lights curves have been transformed into the 
photometric bands UBVRIJK, using bolometric corrections 
and temperature-color relations provided by Castelli, Gratton 
& Kurucz (1997a,b). In this way, light-curve amplitudes and 
mean absolute magnitudes, either intensity-weighted (M) or 
magnitude- weighted (M), are derived in various photometric 
bands. After having tested that intensity- weighted magnitudes 
approximate better the static values than magnitude- weighted 
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M = 0.65Mo 
LogL/Lo=1.91 
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Phase 



Fig. 1. Theoretical bolometric light curves for a subsample (Z=0.001 and M=0.65 M Q ) of fundamental models. The 
luminosity levels are labeled. 



ones, similarly to what found for RR Lyrae (see e.g. Marconi 
et al. 2003) and Classical and Anomalous Cepheids (Caputo, 
Marconi, Ripepi 1999, Marconi et al. 2004), in the following we 
give the predicted relations with ( M) also because empirical 
investigations usually provide this type of mean value^J. 

In Figure [3] we show the behaviour of visual amplitudes 
as a function of the period for different values of metallicity, 
mass and luminosity. As for RR Lyrae models, we find that for 
fixed period the amplitude increases as the luminosity increases 
and as the mass decreases (middle and bottom left panels, re- 
spectively), while it remains quite constant in the considered 
interval of metallicity, at fixed mass and luminosity (top left 
panel). We notice that part of the observed behaviour in the 
period-amplitude diagram is due to the dependence of period 
on mass and luminosity, as simply derived even from a sim- 
ple linear adiabatic approach. The change in the pulsation in 
amplitude is mainly related to the distance from the FBE, as 
shown in the two right panels of the same figure. The deviation 
from linearity of the highest luminosity level models (top panel 
of Figure [3j is related to the complex coupling between pulsa- 
tion and convection for these low density and cool structures, 
(see also Bono et al. 1997c for a similar behaviour in high lu- 

1 Similar relations involving magnitude-weighted mean val- 
ues are available upon request to the authors 



minosity RR Lyrae stars) . However in the range of linearity we 
obtain: 

log P = -0.033(±0.034) - 1.15 log - 0.475 log < M v > 
- O.W5A v (3) 

We notice that this dependence of the period on amplitude 
is, within the errors, in agreement with those found in Di 
Criscienzo et al. 2004. 

4. Same relevant relations 

The linear regression through the data reported in Table 1 
provides the following analytical relations: 

logT e (FBE) = 3.912(±0.005) - 0.035 log + 0.048 log -^-(4) 

Lq Mq 

\ogT e {FRE) = 3.925(±0.005) - 0.075 log — + 0.118 log -£f-(5) 

Lq Mq 

with a rms of 0.005 , valid on the adopted range of metallicity 
and stellar mass. The location of these edges in the My-logP 
plane is presented in Figure [3] and the corresponding analytical 
relations are: 
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Fig. 2. The same as Fig. 2 but for theoretical radial velocity curves. 



< Mv(FBE) > = -0.36(±0.04) - 1.55 log - 2.37 log P 

Mq 



0.05(±0.03)logZ 
0.14(±0.03) - 2.21 
+ 0.05(±0.02)logZ 



<M V (FRE)> = 0.14(±0.03) -2.251og — -2.17 log P 

Mq 



(6) 

y 
(7) 



From the same Figure, we note that the BL Her instability 
strip is the extension at higher luminosities of the fundamental 
instability region for RR Lyrae stars, whereas it does not over- 
lap the one for Anomalous Cepheids (ACs) since the AC FRE 
is hotter than the BL Her FBE (see also Caputo et al. 2004). 
This result, also related to the significant mass difference be- 
tween the two classes of pulsators, is fully consistent with the 
empirical evidence that ACs are brighter than BL Her stars, 
for a given period. 

The natural outcome of the period relation (Eq. 1) in the 
observational plane is the Period-Magnitude-Color relation in 
which the pulsation period is correlated with the pulsator ab- 
solute magnitude and color, for a given mass. As a matter of 
example, linear interpolations through the results give: 



1.10 log 



M 



(My) 



Mq 

-2.08(±0.02) 
M 



(8) 



ilogP-4.13[(M v ) - (M K )} 



2.14 log 



M 



o 



(9) 



with no dependence on metallicity. According to these re- 
lations, for a sample of variables at the same distance and 
with the same reddening, e.g., for variables in a given globu- 
lar cluster or galaxy with no differential reddening, one could 
estimate the mass range spanned by the variables from the ob- 
served luminosity range, once periods and intrinsic colors are 
firmly known. However, these relations have the disadvantage 
of depending on reddening uncertainties and for this reason, 
as discussed in several earlier papers (e.g., Madore 1982, Di 
Criscienzo et al. 2004, Udalski et al.1999), Period- Wesenheit 
relations are widely used, with the Wesenheit functions being 
by definition reddening insensitive. The full set of the predicted 
Period- Wesenheit and Period-Magnitude relations, including 
the evolutionary properties of the BL Her stars, are reported 
in Di Criscienzo et al. (2006). Figure [5] shows the behavior of 
the predicted radii for fundamental pulsators, obtained by av- 
eraging the theoretical radius curve. The solid line is the linear 
regression through the entire set of models (PR relation): 



(M v ) = -1.22(±0.05) -2.60 log P-3.03[(Ms) - (M v )] 



log R = 0.87(±0.01) + 0.529(±0.006) log P 



(10) 
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Fig. 3. (Right panels)Visual amplitudes versus periods for selected models varying the metallicity at fixed mass 
and luminosity (upper panel), varying the luminosity at fixed mass and metallicity (Z=0.001, intermediate panel) 
and varying the mass at fixed luminosity and metallicity (Z=0.001, lower panel). (Left panels) The same but visual 
amplitudes are now plotted versus A log(P)=log (P)- log(PFOB£;) 



where R is in solar units. This relation is in agreement with 
the empirical PR relation obtained by Burki & Meylan, 1986. 

In the same Figure we have also reported the theoretical 
PR relation recently obtained by Marconi et al. (2005) for RR 
Lyrae. This very good agreement between the theoretical PR 
relations obtained for RR Lyrae and BL Her stars supports 
earlier suggestions by Burki & Meylan (1986) concerning the 
similarity of the PR relations for these two classes of Pop. II 
variables, and is consistent with the results by Caputo et al. 
(2004) concerning the behaviour of RR Lyrae and BL Her in 
the Period-Magnitude diagram. 

5. Conclusion and final remarks 

This work is part of a larger project that has the scope to 
study, from the theoretical point of view, Population II vari- 
ables in Globular Clusters and dwarf galaxies. In the last 
few years we have provided the results for RR Lyrae stars 
(Marconi et al 2003, Di Criscienzo et al.2004) and Anomalous 
Cepheids (Marconi, Fiorentino & Caputo, 2004, Caputo et al. 
2004, Fiorentino et al. 2006). The present paper, together with 
the investigation by Di Criscienzo et al. (2006), deals with a 



further class of radial pulsators observed in similar metal poor 
stellar fields: Population II Cepheids. In particular, we present 
the new non-linear and convective pulsational models for short 
period variables, generally named BL Her stars. 

We have investigated the instability strip and shown the 
dependence of the pulsation region on mass and luminosity, 
as well as the effect of varying the the mixing lenght parame- 
ter. We have presented a detailed atlas of light/radial- velocity 
curves. All these curves show a large variety of shapes, which 
is perhaps the most striking feature of BL Her stars as already 
suggested by Moskalik & Buchler (1994) on the basis of non- 
linear but radiative models. 

Theoretical light curves have been transformed in the 
UBVRIJK photometric bands, using the bolometric correc- 
tions and color-temperature relations provided by Castelli, 
Gratton & Kurucz (1997a,b). On this basis, mean magnitudes 
and colors have been evaluated and some relations both for the 
boundaries of the instability strip and for the dependence of 
the absolute magnitude on period, mass and color are derived. 

Finally, we have calculated the period-radius relation and 
found that it is in excellent agreement with the empirical re- 
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Fig. 4. Dependence of the boundaries of the pulsation re- 
gion on luminosity in the My-logP plane for the labeled 
values of mass and metallicity. In the left panel we also re- 
port the AC instability strip for a mass of 1.3M© (Marconi 
et al. 2004) while in the right one we show the predicted 
RR Lyrae instability strip for Z=0.001 as obtained by syn- 
thetic simulations. 
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Fig. 5. Period-Radius relation for our models. The solid 
line is the linear regression throught the entire set of BL 
Her models (filled circle) which is exactly the same found 
by Burki & Meylan, 1996 (BM96), while the dashed one 
is the PR relation obtained by Marconi et al. 2005 (M05) 
for RR Lyrae stars. 



lation by Burki & Meylan and consistent with the one holding 
for shorter periods for RR Lyrae stars. 
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